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Abstract. — Optimal transport enables one to construct a metric on 
the set of (sufRciently small at infinity) probability measures on any (not 
too wild) metric space X. It is called its Wasserstein space 'W2{X) (note 
that there are variant of this space depending on an exponent, but here 
we use only the most common, namely the quadratic cost). 

In this paper we investigate the geometry of y^(X) when X has glob- 
ally non-positive sectional curvature (what we call a Hadamard space). 
Although it is known that -except in the case of the line- 'W2(X) is 
not non-positively curved, our results show that #2(^) have large-scale 
properties reminiscent of that of X. 

First, we are able to define a geodesic boundary in a very similar way 
than one does for X itself, notably thanks to an asymptotic formula that 
gives good control on the distance between two geodesic rays. 

This understanding of the geodesic rays enables us to prove that if X is 
a visibility space (which for example includes all CAT(— 1) spaces, in par- 
ticular all simply connected manifolds with uniformly negative sectional 
curvature), then W2{X) docs not contain any isometrically embedded 
Euclidean plane: one says that it has rank 1. 

We also show that in many cases (including trees and manifolds), if 
X is negatively curved then the isomctry group of W2{X) is as trivial as 
it can be. This contrasts with previous results in the Euclidean case. 

While optimal transport has commonly been studied under lower cur- 
vature bounds (or to define such bounds), we hope that this paper will 
show that it is of interest in a non-positively curved setting. Notably, an 
uncommon feature of our work is that we allow X to be branching, and 
some emphasis is indeed given to trees. 
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1. Introduction 

The goal of this paper is to contribute to the understanding of the 
geometry of Wasserstein spaces. Given a metric space X, the theory 
of optimal transport (with quadratic cost) gives birth to a new metric 
space, made of probability measures on X, often called its Wasserstein 
space and denoted here by W2{X) (precise definitions are recalled in the 
first part of this paper). One can use this theory to study X, for example 
by defining lower Ricci curvature bounds as in the celebrated works of 
Lott-Villani |LV09] and Sturm |Stu06] . Conversely, here we assume 
some understanding of X and try to use it to study geometric properties 
of W2{X). A similar philosophy underlines the works of Lott in |Lot08] 
and Takatsu and Yokota in [TY09j . 

In a previous paper [KloOSj , the second author studied the case when 
X is a Euclidean space. Here we are interested in the far broader class of 
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Hadamard spaces which are roughly the globally non-positively curved 
spaces. The first part of the paper is devoted to recalls on optimal trans- 
port and Hadamard spaces; in particular the precise hypotheses under 
which we shall work are given there (Definition l2.lt see also examples I2.2p . 
Let us stress that we allow X to be branching; trees, product involving 
trees, some buildings are in particular treated in the same framework 
than, for example, symmetric spaces of non-compact type. 

While non-negative curvature is an assumption that is inherited by 
Wasserstein spaces, it is well-known that non-positive curvature is not 
(an argument is recalled in Section 13. 2p . We shall however show that 
some features of Hadamard spaces still hold in their Wasserstein spaces. 
Let us now describe the main results of the article. 

A Hadamard space admits a well-known geometric compactification, 
obtained by adding a boundary at infinity made of asymptote classes of 
geodesic rays. In part HIl we study the geodesic rays of ^(X). Using 
a displacement interpolation procedure (Proposition 14.21) . we associate 
to each ray its asymptotic measure which lies in a subset j3^i{cdX) of 
probability measures on the cone cdX over the boundary of X (Defini- 
tion 15. ip . It encodes the asymptotic distribution of the direction and 
speed of a measure running along the ray. Our first main result is the 
asymptotic formula (Theorem 15.21) which enables one to compute the 
asymptotic behavior of the distance between two rays in terms of the 
Wasserstein distance of the asymptotic measures, with respect to the an- 
gular cone distance on cdX. This asymptotic distance is either bounded 
or asymptotically linear, so that the boundary dW^iX) of the Wasser- 
stein space, defined as the set of asymptote classes of unit geodesic rays, 
inherits an angular metric, just like X does. A striking consequence of the 
asymptotic formula concerns the rank of ^(X), and partially answers a 
question raised in the previous paper cited above. 

Theorem 1.1. — If X is a visibility space (e.g. if it has curvature 
hounded from above by a negative constant), then it is not possible to 
embed the Euclidean plane isometrically in ^(X). 

In other words, when X is strongly negatively curved -which implies 
that its has rank 1-, then although #^(X) is not negatively curved it has 
rank 1 too. Note that our large-scale method in fact implies more general 
non-embedding results, see Proposition 16.31 and the discussion below. 

It is easy to see that any element of !^i{cdX) is the asymptotic mea- 
sure of some ray in #^(X). However, asymptotic measures of complete 
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geodesic must be concentrated on the set of unit geodesies (Proposition 
16. 2p and this is an easy but important step in the proof of Theorem ll.il A 
complete geodesic in /^(X) defines two asymptotic measures, also called 
its ends: one for t — >■ — oo, the other for t — )■ +00. One can therefore 
ask which pairs of boundary measures are the ends of some complete 
geodesic. We give in Section [7] a necessary condition that follows easily 
from the asymptotic formula, and a complete answer when X is a tree 
(Theorem I7.10p . In fact we give two characterizations of realizable pairs 
of asymptotic measure, one of them involving a optimal transportation 
problem related to the Gromov product. This criterion might gener- 
alize to arbitrary Gromov hyperbolic spaces, but the case of a general 
Hadamard space may be challenging. 

Part mil is devoted to the definition of a so-called cone topology on 
dW2{X) and W2{X) = W2{X) U d W2{X) , see Proposition EH Note that 
the angular metric alluded to above, however useful and meaningful, 
does not define a satisfactory topology (just as in dX, where the angular 
metric is usually not separable and can even be discrete). The point is 
that many monotony properties used in the case of Hadamard spaces 
hold when one restricts to angles based at a Dirac mass. This enables 
us to carry out the construction of this topology despite the presence of 
positive curvature. The main result of this part is the following, restated 
as Theorem 19.21 

Theorem 1.2. — The asymptotic measure map defined from dW2{X) 
to ^i{cdX) is a homeomorphism. 

Note that the two natural topologies on dW2{X), namely the cone 
topology and the quotient topology of the topology of uniform conver- 
gence on compact sets, coincide. The set J^i{cdX) is simply endowed 
with the weak topology (where the topology on cdX is induced by the 
cone topology of dX). 

The possibility to identify (9#^(X) to ^i{cdX) should be thought 
of as an interversion result, similar to displacement interpolation. The 
latter says that "a geodesic in the set of measures is a measure on the set 
of geodesies", while the former can be roughly restated as "a boundary 
point of the set of measures is a measure on the (cone over the) set 
of boundary points". Note that #2(X) is not compact; this is quite 
inevitable since W2{X) is not locally compact. 
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Last, Part IIVI adresses the isometry group of #2(-^)- In |Klo08] it 
was proved that the isometry group of y^l^^") is strictly larger than that 
of itself, the case n = 1 being the most striking: some isometries 
of #2(1^) are exotic in the sense that they do not preserve the shape of 
measures. This property seems pretty uncommon, and we shall show the 
following rigidity result. 

Theorem 1.3. — If X is a negatively curved Hadamard manifold or a 
tree not reduced to a line and without leaf, then W2{X) is isometrically 
rigid in the sense that its only isometries those induced by the isometries 
of X itself. 

Note that in the case of trees, the assumption that X has no leaf could 
probably be released but comes out naturally in the proof. We are in 
fact able to treat more space than enclosed in the previous statement, 
see Theorem 110.11 In particular the Ip^q buildings have isometrically rigid 
Wasserstein spaces. 

In the process of proving TLSl we also get the following result (Corollary 
I1U.5I) that seems interesting by itself. 

Theorem 1.4- — Let Y, Z be geodesic Polish spaces and assume that Y 
is geodesically complete and locally compact. Then W2{Y) is isometric to 
W2{Z) if and only ifY is isometric to Z. 

Note that Theorem 11.31 involves the inversion of some kind of Radon 
transform, that seems to be new (for trees it is in particular different 
from the horocycle Radon transform) and could be of interest for other 
problems. 



PART I 
RECALLS AND NOTATIONS 



As its title indicates, this part contains nothing new. We chose to 
give quite a lot of recalls, so that the reader familiar with non-positively 
curved spaces can get a crash-course on Wasserstein spaces, and the 
reader familiar with optimal transport can be introduced to Hadamard 
spaces. 
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2. Hadamard spaces 

We start with some definitions and properties of metric spaces, espe- 
cially tlie non-positively curved ones. Most properties of Hadamard space 
tliat we need are proved in [Bal95] . Another more extensive reference is 
|BH99] . 

2.1. Geodesies. — Let us first fix some conventions for any metric 
space Y (this letter shall be used to design arbitrary spaces, while X 
shall be reserved to the (Hadamard) space under study). 

Given a continuous curve c : I ^ Y, one defines its length by 

k 

i{c) = sup y^ci(Q,_,,Qj g [0, +00] 

where the supremum is over sequences of arbitrary size of tj G /. The 
curve is said to be rectifiable if it has finite length. 

A geodesic in F is a curve 7 : / — )■ y defined on some interval /, such 
that there is a constant v that makes the following hold for all times t, t': 

d{jt,lt') = v\t - t'\. 

In particular, all geodesies are assumed to be globally minimizing and to 
have constant, non necessarily unitary speed. Of course, any geodesic 
defined over a segment is rectifiable, and has length equal to the distance 
between its endpoints. When w = we say that the geodesic is constant 
and it will be necessary to consider this case. We denote by ^^'^'(F) the 
set of geodesies defined on the interval [T, T'] . A geodesic ray (or ray, or 
complete ray) is a geodesic defined on the interval [0,+oo). A complete 
geodesic is a geodesic defined on M. The set of rays is denoted by M{Y), 
the set of unit speed rays by H^iiY) and the set of non-constant rays by 
MyQ{Y). We shall also denote by $f'*(F) the set of complete geodesies, 
and by '^f{Y) the set of unit-speed complete geodesies. 

2.2. Non-positive curvature. — A metric space is geodesic if any 
pair of points can be linked by a geodesic. A triangle in a geodesic 
space Y is the data of three points (x, y, z) together with three geodesies 
parametrized on [0, 1] linking x to y, y to z and z to x. Given a tri- 
angle, one defines its comparison triangle (x, y, z) as any triangle of the 
Euclidean plane that has the same side lengths: d{x,y) = d{x,y), 
d{x,z) = d{x,z) and d{y,z) = d{y,z). The comparison triangle is de- 
fined up to congruence. 
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A triangle with vertices {x, y, z) is said to satisfy the CAT(O) inequahty 
along its [yz\ side (parametrized by a geodesic 7 G ?f '^(F)) if for all 
t G [0, 1], the following inequality holds: 

(1) d{x,-i{t)) ^d{x,{l-t)y + tz) 

see figure [H A geodesic space is said to be locally CAT(O) if every point 
admits a neighborhood where all triangles satisfy the CAT(O) inequality 
(along all there sides). When F is a Riemannian manifold, this is equiva- 
lent to ask that Y has non-positive sectional curvature. A geodesic space 
is said to be globally CAT(O) if all its triangles satisfy the CAT(O) in- 
equality. Globally CAT(O) is equivalent to simply connected plus locally 
CAT(O). We shall simply say CAT(O) for "globally CAT(O)", but this is 
not a universal convention. 




Figure 1. The CAT(O) inequality: the dashed segment is 
shorter in the original triangle on the left than in the com- 
parison triangle on the right. 



2.2.1. Angles. — The CAT(O) condition can be translated in terms of 
angles as follows. Given any geodesic triangle, choose any of its vertices, 
say a;, and assume that the sides containing x are parametrized by two 
geodesies a, 7 G ^°'^(F). If Y is CAT(O), then the Euchdean angle ^^^^xat 
at X is a nondecreasing function of s and t. 

One then defines in Y the angle Z-yixo-i at x as the limit, when s and 
t go to zero, of Z^^^^^. As a consequence, one gets that for any geodesic 
triangle with vertices (x, z) in a CAT(O) space, angle and comparison 
angle satisfy 

2.2.2. Distance convexity. — In a CAT(O) space, given two geodesies 7 
and /3, the distance function t 1 — y d{'-ft, A) is convex. This important 
property shall be kept in mind since it will be used very often in the 
sequel. 
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2.3. Hadamard spaces: definition and examples. — We can now 

introduce the class of spaces we are interested in. 

Definition 2.1. — A metric space is a Hadamard space if it is: 

— Polish (i.e. complete and separable), 

— locally compact, 

— geodesic, 

— CAT(O), implying that it is simply connected. 

In all what follows, we consider a Hadamard space X. The Hadamard 
assumption may not always be made explicit, but the use of the letter X 
shall always implicitely imply it. 

Example 2.2. — There are many important examples of Hadamard 
spaces. Let us give some of them: 

— the Euchdean space W^, 

— the real hyperbolic space MH", 

— the other hyperbolic spaces CH", HH", OH^, 

— more generally the symmetric spaces of non-compact type, like the 
quotient SL(n; R)/SO(n) endowed with the metric induced by the Killing 
form of SL(n; M), 

— more generally any simply connected Riemannian manifold whose 
sectional curvature is non-positive, 

— trees, 

— any product of Hadamard spaces, 

— some buildings, like product of trees having unit edges and no leaf 
or Ipq buildings (see |Bou97|. |BP99]), 

— the gluing of any two Hadamard spaces along isometric, convex sub- 
sets; for example any Hadamard space with an additional geodesic ray 
glued at some point, or three hyperbolic half-planes glued along their 
limiting geodesic, etc. 

2.4. Geodesic boundary. — The construction of the geodesic bound- 
ary that we will shortly describe seems to date back to jEQ73] . but note 
that |Bus55] is at the origin of many related ideas. 

Two rays of X are asymptotic if they stay at bounded distance when 
t — 7- +00, and this relation is denoted by ~. The asymptote class of a ray 
7 is often denoted by 7(00) or 700, and is called the endpoint or boundary 
point of 7. 
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The geodesic (or Hadamard) boundary of X is defined as the set 



Using the convexity of distance along geodesies, one can for example 
prove that, given points x G X and C G dX, there is a unique unit ray 
starting at x and ending at C,. 

The union X = XVJ dX can be endowed with its so-called cone topol- 
ogy, which makes X and dX compact. Without entering into the details, 
let us say that this topology induces the original topology on X, and that 
given a base point Xq a basic neighborhood of a point C = 7(00) G dX 
(where 7 starts at xq) is the union, over all rays a such that d{at^ 7t) < e 
for all t < R, oi a{[R, +00]) (see figure [2]). 



Figure 2. A basic neighborhood of a point CdX, in the cone topology. 

Consistently whith the cone topology, all previously defined sets of 
geodesies, as well as the larger sets C{I,X) of continuous curves defined 
on an interval / with values in X, are endowed with the topology of 
uniform convergence on compact sets. With this topology, since X is 
Hadamard, a geodesic segment is uniquely and continuously defined by 
its endpoints and a geodesic ray is uniquely and continuously defined 
by its starting point, its endpoint in the boundary and its speed. As a 
consequence, there are natural homeomorphisms 



where cdX is the cone over dX, that is the quotient of dX x [0, +00) 
by the relation (x, 0) ~ (y, 0) for all x,y E dX. We usually use the 
same notation (x, s) for a couple and its equivalence class under this 
relation; here s shall be considered as a speed. In particular one has 



Note that in view of our assumptions on X, all these spaces are locally 
compact Polish topological spaces (that is, metrizable by a metric that 



dX = ^i{X)/ 




'^^'^'{X) ~ X^, ^{X) ~ X X cdX 



{M{X)/ ~) ~ cdX. 
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is separable and complete). It ensures that all finite measures on these 
spaces are Radon. 

2.5. Possible additional assumptions. — At some points, we shall 
make explicitely additional hypotheses on X . One says that a space is: 

— geodesically complete if every geodesic can be extended to a complete 
geodesic, 

— non-branching if two geodesies that coincide on an open set of pa- 
rameters coincide everywhere on their common definition interval, 

— CAT(k) fo some k < if its triangles satisfy the ([T]) inequality 
when the comparison triangle is taken in MH^, the hyperbolic plane of 
curvature k, instead of M^. 

— a visibility space if for all pairs of distinct points a,uj E dX there is 
a complete geodesic 7 such that •y^—oo) = a and 7(+oo) = u. 

Note that for all k < 0, the condition CAT(k) implies that X is a visibility 
space. 

Geodesic completeness is quite mild (avoiding examples as trees with 
leafs) , while the other are strong assumptions (for example non-branching 
rules out trees and visibility rules out products). Another possible as- 
sumption, that we shall not use directly, is for X to have rank one, 
meaning that it admits no isometric embedding of the Euclidean plane. 
It is a weaker condition than visibility. More generally the rank of a 
space Y is the maximal dimension of an embedded Eulidean space, and 
it has been proved a very important invariant in the study of symmetric 
spaces. For example the hyperbolic spaces MH*^, CH", HH" and OH^ are 
the only rank one symmetric spaces of non-compact type. 

2.6. Asymptotic distance. — Given two rays 7, a, one defines their 
asymptotic distance by 



This limit always exists because of the convexity of the distance function 
along geodesies. Moreover d^o defines a metric on cdX and, by restriction, 
on dX (in particular two rays whose distance grows sub-linearly must 
be asymptotic). It can be proved that d^o is the cone metric over dX 
endowed with the angular metric. Namely, for any (^,s), (^', t) G cdX, 




t— s>+oo 



t 





GEOMETRIC WASSERSTEIN SPACES 



11 



where Z(,^, C,') = sup^^x ^^(C, is the supremum of angles between rays 
issuing from x and asympotic to and C,' respectively (we refer to [BHQQj 
Section II. 9] for more details and proof). 

It is most important to keep in mind that the metric d^o does not 
induce the cone topology, but a much more finer topology. The most 
extreme case is that of visibility spaces, where doo{'~f, a) is if 7 ~ a and 
the sum of the speeds of 7 and cr otherwise: the topology induced on 
dX is discrete. However, the function d^o is lower semi-continuous with 
respect to the cone topology |BH99t Proposition II. 9. 5], so that it is a 
mesurable function (see the next paragraph for a reference). 

In higher rank spaces, it can be useful to turn d^o into a length metric, 
called the Tits metric, but we shall not use it so we refer the interested 
reader to the books cited above. Let us just note that doo resembles in 
some aspects the chordal metric on a sphere. In particular, it is naturally 
isometric to this metric when X is a Euclidean space. 



In this section, we recall the definition of Wassertein space and some 
of its main properties. For more details, we refer to the books |Vil03] 
and |Vil09] . 

3.1. Optimal transport and Wasserstein metric. — Let us start 
with the concept of optimal transport which is the theory aimed at study- 
ing the Monge-Kantorovich problem. 

3.1.1. Optimal mass transport problem. — Standard data for this prob- 
lem are the following. We are given a Polish metric space (F, d), a lower 
semicontinuous and nonnegative function c : F x F — called the cost 
function and two Borel probability measures /i, v defined on F. A trans- 
port plan between /i and is a measure on y x y whose marginals are /x 
and u. One should think of a transport plan as a specification of how the 
mass in Y , distributed according to /i, is moved so as to be distributed 
according to v. We denote by r(yU, v) the set of transport plans which is 
never empty (it contains fi^u) and most of the time not reduced to one 
element. The Monge-Kantorovich problem is now 



3. The Wassertein space 




12 



JEROME BERTRAND & BENOIT KLOECKNER 



where a minimizer is called an optimal transport plan. The set of optimal 
transport plans is denoted by To{^,cr)- 

Let us make a few comments on this problem. First, note that un- 
der these assumptions, the cost function is measurable (see, for instance, 
jVil031 p. 26]). Now, existence of minimizers follows readily from the 
lower semicontinuity of the cost function together with the following com- 
pactness result which will be used throughout this paper. We refer to 
jBil99] for a proof. 

Theorem 3.1 (Prokhorov's Theorem). — Given a Polish space Y , 
a set P C l3^{Y) is totally hounded (that is, has compact closure) for the 
weak topology if and only if it is tight, namely for any e > 0, there exists 
a compact set such that fi(Y \ K^) < e for any e P. 

For example, the set r(/i, v) is always compact. 

We also mention that, compared to the existence problem, the issue 
of the uniqueness of minimizers is considerably harder and requires, in 
general, additional assumptions. To conclude this introduction, we state 
a useful criterium to detect optimal transport plan among other plans, 
named cyclical monotonicity. 

Definition 3.2 (cyclical monotonocity). — Given a cost function 
c : Y X Y — y ]R+, a set r C F X y is called c-cyclically monotone 
if for any finite family of pairs (xi, yi), - ■ ■ , {xm, ym) in the following 
inequality holds 

m m 

(3) ^ c(xi, yi+i) > ^ c(xi, yi) 

i=l i=l 

where ym+i = yi- 

In other words, a c-cyclically monotone set does not contain cycles 
of pairs (starting point, ending point) along which a shift in the ending 
points would reduce the total cost. 

Theorem 3.3. — Let (Y, d) he a Polish space and c:Y xY — )■ R"*" he 
a lower semi- continuous cost function. Then, a transport plan is optimal 
relatively to c if and only if it is concentrated on a c-cyclically monotone 
set. 

If c is continuous, this is equivalent to its support heing c-cyclically 
monotone. 
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Under these assumptions, this result is due to Schachermayer and Te- 
ichmann |ST09j : see also |Vil09] for a proof. 

3.1.2. Wassertein space. — Wasserstein spaces arise in a particular vari- 
ant of the setting above. 

Definition 3-4 (Wassertein space). — Given a Polish metric space 
Y, its (quadratic) Wassertein space W2(Y) is the set of Borel probability 
measures fi on Y with finite second moment, that is such that 



J d{xo, xY fi{dx) < +00 for some, hence all Xq G Y, 
endowed with the Wasserstein metric defined by 



From now on, the cost c will therefore he c = d^. 
The fact that W is indeed a metric follows from the so-called "gluing 
lemma" which enables one to propagate the triangular inequality, see e.g. 



Remark 3.5. — In this paper, we will also consider the Wassertein 
space over the cone cdX relative to the cost d^. Since d^ is lower semi- 
continuous, this is indeed a metric space. Note that we shall denote by 
Woo the Wasserstein metric derived from doo- 

The Wasserstein space has several nice properties: it is Polish; it is 
compact as soon as Y is, in which case the Wasserstein metric metrizes 
the weak topology; but if Y is not compact, then W2{Y) is not even 
locally compact and the Wasserstein metric induces a topology stronger 
than the weak one (more precisely, convergence in Wasserstein distance is 
equivalent to weak convergence plus convergence of the second moment). 
A very important property is that /^(F) is geodesic as soon as Y is; let 
us give some details. 

3.2. Displacement interpolation. — The proof of what we explain 
now can be found for example in chapter 7 of |Vil09] . see in particular 
corollary 7.22 and Theorem 7.30. Note that the concept of displacement 
interpolation has been introduced by McCann in |McC97] . We write 
this section in the case of X, but most of it stays true for any Polish 
geodesic space. 




|Vil09] . 
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3.2.1. Dynamical transport plans. — Define a dynamical transport plan 
between two measures /io,/ii G W2{X) as a probability measure /i on 
C{[to, ti]] X) such that for i = 0, 1 the law at time ti of a random curve 
drawn with law /i is //j. In other words we ask Cj^^/i = /ij where Ct is the 
map C{[tQ,ti]; X) — >■ X defined by 6^(7) = ■jf 

The cosi of /i is then 

where £(7) is the length of the curve 7 (possibly +00). A dynamical trans- 
port plan is optimal if it minimizes the cost over all dynamical transport 
plans, and it is known that a dynamical transport plan exists. Morevover, 
if /i is an optimal dynamical transport plan then: 

(1) the law (ej(,,efj#/i of the couple (7jo,7ti) where 7 is a random 
curve drawn with law n, is an optimal transport plan between /iq and /ii, 

(2) /i-almost all 7 G C{[tQ,ti]; X) are geodesies. 

Conversely, if 11 is a (non-dynamical) optimal transport plan, then one 
can construct for any < ^1 an optimal dynamical transport plan by 
the following construction. Let F : X"^ — )■ C{[tQ,ti]; X) be the map that 
sends a couple (x, y) of points to the unique geodesic parametrized on 
[^0,^1] that starts at x and ends at y. Then /i = F^U is an optimal 
dynamical transport plan, whose associated optimal transport plan is 
obviously 11. 

Given a dynamical plan /i G ^(^'^(X)) where / is an arbitrary inter- 
val, /i is said to be c-cyclically monotone if for any s,t E I, the support 
of the plan [cg, et)^fi is c-cyclically monotone (note that here the cost is 
continuous). As soon as fit have finite second moments, this is equivalent 
to being optimal, but cychcal monotonicity has the advantage of being 
well-defined without any integrability assumption. 

3.2.2. Consequences of optimality. — The main use of optimal dynam- 
ical transport plans is that they define geodesic segments. Indeed, let as 
before be the map 7 1— )■ 7t defined on the set of continuous curves. If 
/i is an optimal dynamical transport plan, consider the law /i^ = Cj^/i 
at time t of a random geodesic drawn with law fi: then {fit)tQ^t^ti is a 
geodesic of y^(X). Displacement interpolation is the converse to this 
principle. 
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Proposition 3.6 (Displacement interpolation) 

Given any geodesic segment (/it)to=£t^ti ^(-^)> there is a probability 
measure fi on such that for all t fit = Ct^fi. 



If /i is a dynamical transport plan on /, for all to)^i ^ ^ define the 
time restriction of /i to [to,^i] as = rtf^^n^l^ where ^(^^(^(7) is the 

restriction of the curve 7 to the interval [to;^i]- 

Let /i be an optimal dynamical transport plan on [0, 1]. For all to, ti G 
[0, 1] the following hold (see |Vil09[ Theorem 7.30]): 



(1) is an optimal dynamical transport plan, 

(2) if X is non-branching and (to,ti) 7^ (0, 1), then is the unique 
(up to parametrization) optimal dynamical transport plan between fito 
and fiti- 

3.2.3. The Wasserstein space is not non-positively curved. — It is well 
known that the non-positive curvature assumption on X is not inherited 
by W2{X) except if X is the real line or a subset of it. Let us use 
displacement interpolation to show it in the most simple way. First, the 
Wasserstein space is contractible thanks to the affine structure on the set 
of probablity measures (the convex combination of fi and u with weights 
t and (1 — t) being tfi+ {1 — t)u). If it were locally CAT(O), it would be 
globally CAT(O) and, in particular, uniquely geodesic. But as soon as 
X is not reduced to a geodesic, there exists four distinct points x, x', z 
such that d{x,y) = d{x,z) and d{x',y) = d{x',z) (see figure [3] for the 
case of a tree). Between the measures /i = ^6x + (where 5's are 
Dirac masses) and z/ = ^6y + ^6z, all transport plans are optimal (and 
there are several of them, in fact an uncountable number, for example 
parametrized by the quantity of mass taken from x to y). As we have 
seen above, each optimal transport plan defines a geodesic in #2(^) from 
fi to u. Therefore #2(^) is very far from being uniquely geodesic, and is 
in particular not non-positively curved. 



16 



JEROME BERTRAND & BENOIT KLOECKNER 




Figure 3. Between a measure supported on x and x' and a 
measure supported on y and z, there are many optimal trans- 
port plans and therefore many geodesies. 

PART II 
GEODESIC RAYS 

In this part, we study the geodesic rays of W2{X), which shall in par- 
ticular be used to define its geodesic boundary, to be endowed with a 
topology in the next part. 

The geodesic boundary of W2{X) is defined as the set of asymptote 
classes of unitary geodesic rays: 

dW2{X)=Mi{W2{X))/ ^ . 

The core result of this part is the asymptotic formula (Theorem 15.21) 
which gives precise indication on the speed at which the distance between 
two rays goes to infinity. In particular, we shall show that this distance 
is either bounded or linear, so that there is a well-defined asymptotic 
distance in The asymptotic formula shows that it can be iden- 

tified with the Wasserstein metric derived from the asymptotic distance 
of X. 

To do this, we need to extend the displacement interpolation principle 
to rays (Section Hj), so that we can define the asymptotic measure of a 
ray that yields the identification alluded to above (Section [5]). Then we 
use the asymptotic formula to study the rank of #2(^) (Section E]), and 
last we discuss the possibility to prescribe the asymptotic behavior of a 
complete geodesic for t — t- +oo and t — t- — oo at the same time (Section 

ED- 
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4. Displacement interpolation for rays 

To start the description of we need displacement interpolation 

for rays. There is not much, but a little work needed to extend from the 
case of geodesic segment, the crucial point being to handle the branching 
case. Note that the case of complete geodesic is not different than that 
of the rays. 

Proposition (Displacement interpolation for rays in the 
non-branching case) 

If X is n on- branching, any geodesic ray (fit) o/#2(-^) admits a unique 
displacement interpolation, that is a probability measure fi on ^%{X) such 
that Ht is the law of the time t of a random ray having law fi (in other 
words, such that fit = {et)#fi). 

Note that since (/ij) is a geodesic, the time restriction of fi to any 
segment is an optimal dynamical transport plan. 

Proof. — Since X is non-branching and (/Xt) is defined for all positive 
times, we know that there is a unique optimal dynamical transport plan 
/i(T) from /io to ht, and that fl{T) = jl{T')^'^ whenever T < T'. The fact 
that X is non-branching also implies that for any two points x,y & X 
there is a unique maximal geodesic ray F{x, y) starting at x and passing 
at y at time 1. Moreover this ray depends continuously and therefore 
measurably on {x,y). It follows that all /i(T) are uniquely defined by 
/i(l) and that the probability measure /i = F^{Il), where 11 is the optimal 
transport plan associated to has the requiered property. □ 

When X is branching, the previous proof fails for two reasons. The 
first one is that we cannot determine fl{T) and define /i from /i(l) alone; 
although a theorem of Prokhorov will do the trick. The second problem 
is that there may exist several optimal dynamical transport plans corre- 
sponding to the same geodesic; but the set of these transports is always 
compact and a diagonal process will solve the problem. However, we lose 
unicity in the process and it would be interesting to single out one of the 
dynamical transport plan obtained. 

Proposition 4-2 (Displacement interpolation for rays) 

Any geodesic ray of W2{X) admits a displacement interpolation. 

Proof — Let {fit)t^o be a geodesic ray in W2{X) and for all T > 0, let 
M{T) be the set of all optimal dynamical transport plans parametrized 



18 



JEROME BERTRAND & BENOIT KLOECKNER 



on [0,T] that induce the geodesic segment (yUt)o^t^T- It is a compact 
set according to |Vil09] . Corollary 7.22. For all T G N, choose /i(T) in 
M(T) and for all integer ^ T' ^ T define fi{T'\T) as the restriction 

By a diagonal process, one can extract an increasing sequence of 
integers such that for all T' e N, fi{T'\Tk) has a limit /i(T') e M(T') 
when k — )■ +oo. Since for all T' < T the restriction maps 

are continuous, we get that p^'^' (fL(T)) = fi{T'). We therefore have a 
projective system of measures; the projection map 

commutes with the p'^''^' thus, according to a theorem of Prokhorov (see 
|Sch70] ). if we prove tightness, i.e. that for all e > there is a compact 
K C ^{X) such that for all T, jl{T){p^ K) ^ 1 — e, then we can conclude 
that there is a unique measure ^ on ^(X) such that p^/x = /i(T) for all 
T. This measure will have the requiered property since Ji{T) G M(T). 

Fix any e > 0. Let Kq, Ki be compact subsets of X such that fii{Ki) ^ 
1 — £:/2 for z = 0, 1. Let K be the compact subset of M{X) consisting in 
all geodesic rays starting in Kq and whose time 1 is in Ki. Then for all 
T > 1, ^{T){p^{K)) ^ 1 - e, as needed. □ 

The following result shall make displacement interpolation particularly 
useful. 

Lemma 4-3 (lifting). — Let fi, a be probability measures onM{X) (or 
similarly ^^'^ (X), ...) and denote by Ht = {,^t)#l^ o^nd at = {et)jj,(J their 
time t. 

Any transport plan G r(/ij, ctj) admits a lift, that is a transport plan 
n G r(/i, a) such that lit = (e^, eJ^LI. 

Proof. — Disintegrate /i along fit- there is a family (Cx)xgx of probability 
measures on ^(X), each one supported on the set ejr^(x) of geodesic rays 
passing at x at time t, such that yU = / CxfJ't{dx) in the sense that 

/x(A) = / Cx(^)/^t(c?2;) 

for all measurable A. Similarly, write cr = J ^yCTt{dy) the disintegration 
of a along ctj. 
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Define then 

U{A xB)= f UA)^y{B)Ut{dxdy). 

It is a probability measure on ^(X)^, and for any measurable sets A, B 
in X we have Ii{et^{A), eJ^{B)) = Ut{A x B) because (x {eT\A)) is 1 if 
X G A, otherwise (and similarly for (e^^(i?))). A similar computation 
gives that 11 has marginals fi and a. □ 

Note that we gave the proof for the sake of completness, but one could 
simply apply twice the gluing lemma, after noticing that the projection 
et from ^{X) to X gives deterministic transport plans in r(/i, /xj) and 

The lift of lit need not be unique; the one constructed in the proof is 
very peculiar, and can be called the most independent lift of IIj. It is well 
defined since other disintegration families {C'x)x and {^y)y must coincide 
with {Cx)x and {^y)y for /xt-almost all x and aj-almost all y respectively. 

The lifting lemma shall be used to translate the optimal transport 
problems between fit and at, where these measures move (usually along 
geodesies) to transport problems between the fixed and a, where it 
is the cost that moves. In other words, we have just shown that mini- 
mizing J c{x,y) Ut over the Ut G r{fif,<^t) is the same than minimizing 
/ c(7t, /3t) U{d-fd^) over the U e r(/i, a). 



5. Asymptotic measures 

Let us denote by Coo the map defined by the formula 

. .^ Coo : ^{X) — y cdX 

^ ^ 7 ^ (M,^(7)) 

where 7^ is the unitary reparametrization of 7, [7^] is its asymptote class 
and 5(7) is the speed of 7. It is to be understood that whenever 3(7) = 0, 
[7^] can be taken arbitrarily in dX and this choice does not matter. 

Definition 5.1 (asymptotic measure). — Let {fit)t>o be a geodesic 
ray in y^(X) and /i be a displacement interpolation (so that fit = ^t#f^)- 
We define the asymptotic measure of the ray by 
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We denote by ^i{cdX) the set of probability measures u on cdX such 
that / v^u{dv) = 1. 

In the branching case, the dynamical optimal transport plan is not 
unique in general. Therefore, the asymptotic measure depends a priori 
on the choice of the dynamical optimal transport plan. We will see soon 
that it is not the case. 

Note that the speed of the geodesic (fit) is 

\ 1/2 / \ 1/2 



and we denote it by s{fi). In particular, J^i{cdX) is the set of measures 
that correspond to unit speed geodesies. We shall use that the speed 
function s defined in M{X) is in several times. 

The main result of this section is the following. 

Theorem 5.2 (asymptotic formula). — Consider two geodesic rays 
{l^t)t>o and {at)t>o, let fi and a be any of their displacement interpolations 
and fioo, CToo be the corresponding asymptotic measures. Then (fit) and 
(at) are asymptotic if and only if fioo = <^oo, o,nd we have 

hm = Woo(/ioo,o-oo). 

t— >oo t 

Therefore: as in X itself, the distance between two rays is either 
bounded or of linear growth, and two displacement interpolations of the 
same ray define the same asymptotic measure. 

The role of the asymptotic formula goes far beyond justifying Def- 
inition 15.11 in the branching case: it gives us a very good control on 
geodesic rays of W2{X) on which several of our results rely. To cite one, 
the asymptotic formula is the main ingredient of Theorem 11.11 on the 
rank oiW2{X). 

For every t ^ 0, let dt be the function defined on ^(X) x M{X) 
by dt{pf,(3) = d{'jt,/3t)- We start with an implementation of a classical 
principle. 

Lemma 5.3. — The function dt is in L^(r(/i, cr)) , by which we mean 
that there is a constant C = C{fi,a) such that for all U G r{fi,a), 

Jd^n^c. 

In the following, it will be of primary importance that C does not 
depend on 11. 
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Proof. 



Denoting by x any base point in X we have 




□ 



Proof of the asymptotic formula. — We shall use several time the fol- 
lowing translation of the convexity of the distance function: given any 
geodesic rays 7, (3 of X, the function 



is nondecreasing in t and has limit d^{'y, 

Assume first that fi^o = o"oo, and let us prove that W(/^t, is bounded. 
The lifting lemma gives us a transport plan 11 G r(/i, a) such that for 
all (7, f3) in its support, 700 = /3oo (simply lift the trivial transport (Id x 
Id)#/ioo)- If ft{l,f3) is positive for some t, then doo{l,P) > 0. It follows 
that on suppll, dt ^ do. Therefore: 



which is bounded by the previous lemma. 

Let now Uoo be a transport plan from fi^o to a^o that is optimal with 
respect to d^o. Such a minimizer exists since d^o is non-negative and 
lower-semicontinuous with respect to the cone topology of cdX; note 
that taking an almost minimizer would be sufficient anyway. Denote by 
n a lift of IIoo to r(/i, a); then 



We have 2 ^ (ioo ^ /< ^ /i ^ —do for all t ^ 1, so that t^^dt is bounded 
hj 2 + do and —do. We can thus apply the dominated convergence theo- 
rem, which gives 



/*(7,/3): 



d{-ft, A) - d{-fo, /3o) 



t 
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To prove the other inequahty, we introduce gt := max(0, /j). It is a 
nondecreasing, nonnegative function with doo as hmit, and it satisfies 

t'g^ ^ dl 

Let lit be an optimal transport plan between (/if) and [at) and lit be 
a lift to r(/i,(T), which by Prokhorov Theorem is compact in the weak 
topology. Let {tk)k be an increasing sequence such that 

limW(/itfe,a-fJ = liminf W(/Xt,o-t) 

k t 

and Ilf^ weakly converges to some Iloo- 
For all k' < k, we have 

W'K,afJ _ 1 




Letting k ^ oo and using that the gt are continuous, we obtain 

hm^mf ^ / ^f^, Hoo 

for all k'. But (/f^, ^ 2 and the dominated convergence theorem enables 
us to let k' — oo: 

Immf ^ d^Uoo^ d^{Hoo,(Toc)- 

This ends the proof of the asymptotic formula, and shows that if (/if) 
and (cTf) stay at bounded distance then /ioo = Coo- CH 

The end of this part consists in the investigation of two questions where 
the asymptotic formula finds applications, but from which the sequel is 
independent. The reader interested in the topology of W2{X) can go 
directly to Part IIIII 

6. Complete geodesies and the rank 

In this section we study complete geodesies in /^(X), in particular to 
understand its rank. Recall that the rank of a metric space is the highest 
dimension of a Euclidean space that embeds isometrically in it. 

The main result of this section is the following. 



/ < Hf, 
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Theorem 6.1. — If X is a visibility space, in particular if it is CAT(k) 
with K, < 0, then y^(X) has rank 1. 

We expect more generally that for all Hadamard space X, the rank 
of #^(X) is equal to the rank of X, or at least bounded above by some 
weak rank of X. Theorem 16.11 is a first step in this direction. Note 
that the fact that we ask the embedding to be nothing weaker than an 
isometry and to be global is important: it is easy to prove that as soon 
as X contains a complete geodesic, it admits bi-Lipschitz embedding 
of Euclidean spaces of arbitrary dimension, and isometric embeddings of 
Euclidean balls of arbitrary radius and dimension. Both these facts come 
from an isometric embeeding of open Euclidean half-cone of arbitrary 
dimension into #2(1^)5 see |Klo08] . 

First, let us see that the asymptotic measure of a complete geodesic is 
much more constraint than that of a mere ray. 

Proposition 6.2. — Let {fit)teM. be a complete unit speed geodesic of 
#2(^)7 di^d fi be one of its displacement interpolations. Then fi is con- 
centrated on the set of unit speed geodesies of X . 

Note that here we do not use any assumption on X (besides the exis- 
tence of displacement interpolations). 

Proof — Let 7, f3 be two geodesies in the support of /i and let a = 3(7) 
and b = s{(3). Fix some point x G X. Then we have the equivalents 
d{x, '-ft) ~ at and d{x, Pt) ~ bt when t — )■ ±00. In particular, we get that 



and similarly d'^{(3t,'y-t) ^ (a + b^t'^ + o(t^). We also have d'^ijttl-t) = 
AaH"^ and d'^{[3t,P-t) = 4h'^t^. But, since is the displacement interpola- 
tion of a geodesic, the transport plan 11* it induces from /if to H-t must 
respect the cyclical monotonicity. In particular we have 



From this and letting t — 00 we get 4a^ + 46^ ^ 2a^ + 26^ + 4a6, which 
is only possible when a = b. 

We proved that the speed of geodesies in the support of fi is constant, 
and since their square integrates (with respect to fi) to 1 we get the 



d\-ft,f3_t) ^ {d{-ft,x) + d{x,f3^t)f 
^ {a + b)H^ + o{t^) 



d\jt,l-t) + d\f3t,f3^t) <: d\jt,f3-t) + d^ 



(A,7-<). 



desired conclusion. 



□ 
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As a consequence, the assymptotic measure of a ray that can be ex- 
tended to a complete geodesic hes in the subset ^{dX) of ^^i{cdX) 
(where we identify a space Y with the level F x {1} in its cone). 

Proposition 6.3. — If X is a visibility space, the space ^{dX) en- 
dowed with the metric Woo ( where a ray is identified with its asymptotic 
measure) contains no non-constant rectifiable curve. 

Proof. — First, the visibility assumption implies that whenever 7 and P 
are non asymptotic, unit speed geodesic rays of X, we have (ioo(7, (3) = 2. 
It follows that for all displacement interpolation fj,,a of rays in W2{X), 
assumed to be concentrated on Mi{X), we have 

(5) Woo(/ioo,0-oo) = 2|/ioo - (Tooiy^ 

where | ■ |„ is the total variation norm. 

The proposition now results from the more general following lemma, 
which is well-known at least in the case of Euclidean space. 

Lemma 6.4 (Snowfiaked metrics). — Let{Y,d) be any metric space 
and a < 1 be a positive number. Then {Y, d°') is a metric space not 
containing any non-constant rectifiable curve. 

The fact that d" is a metric comes from the inequality a + 6 ^ (a" + 
^ay/a positive a, 6. 

Let c : / — y be a non-constant curve. Up to restriction and 
reparametrization, we can assume that / = [0, 1] and cq 7^ ci. Take 
any positive integer n; since d is continuous, by the intermediate value 
theorem there are numbers ti = < ^2 < ■ ■ ■ < < 1 such that 
d{ct,_^,ctj = d{co,Ci)/n and d{ct^,Ci) ^ d{co,Ci)/n. Denoting by i the 
length according to the "snowfiaked" metric rf", we get that 



Proof of Theorem \6.1\ — Assume that there is an isometric embedding 
(y9 : — )■ #2(^)- Let be the ray starting at the origin and making an 
angle 9 with some fixed direction. Then, since extends to a complete 
geodesic, so does cp o r^. The displacement interpolation of this ray 
of W2{X) must be concentrated on Mi{X) by Proposition 16. 2^ so that 
/i^ G /^{dX). But (p being isometric, the map 9 — )■ should be an 
isometric embedding from the boundary of (that is, the unit circle 




Since this holds for all n, i{c) = 00 and c is not rectifiable. 



□ 
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endowed with the chordal metric) to (^((9X), Woo)- In particular its 
image would be a non-constant rectifiable curve, in contradiction with 



To prove Theorem I6.1[ one could also simply consider the three Eu- 
clidean geodesies depicted in figure S] and check that their relative dis- 
tances cannot be realized in ^{dX) endowed with Woo- However, our 
method yields easily more general results: to give two simple examples. 
Proposition 16.31 rules out the isometric embedding of Minkowski planes 
(R^ endowed with any norm), and their cones of the form {x^ < ey"^} for 
any e > 0. This contrast with above-mentionned fact that even when X 
is reduced to a line, some Euclidean half cones of arbitrary dimension 
embeds isometrically in y^(X). 



Figure 4. A subset of the plane that cannot appear even 
asymptotically in W2{X). 



7. Pairs of boundary measures joined by a geodesic 

A complete geodesic (/it) in #^(X) defines two rays (for t ^ and 
t ^ 0), and one therefore gets two asymptotic measures, denoted by 
/i_oo and /U.+005 also called the ends of the geodesic. By Proposition 16. 2[ 
these measures are concentrated on dX viewed as a subset of cdX. In 
particular, #2(-^) is already far from being a visibility space. One can 
ask how far precisely: 



Proposition 16.31 



□ 




V2 - V2t 



Question 1. — Which pairs of probability measures on dX are the ends 
of some complete geodesic? 



26 



JEROME BERTRAND & BENOIT KLOECKNER 



Note that we shall need to consider measures /i on the set of unit 
complete geodesies ^"^{X) that satisfy the cyclical monotonicity, but 
such that Ct^/i need not have finite second moment. We still call such 
maps dynamical transport plan and we say that e±oo#/^ are its ends. 
Such a measure defines a complete unit geodesic in if and only if 

€ W2{X) for some, hence all t G M. Note that in this section, we 
only consider unit geodesies even if it is not stated explicitely. 

7.1. A first necessary condition: antipodality. — The asymptotic 
formula shall give us a first necessary condition. 

Definition 7.1 (antipodality). — Two points C;'C ^ "^X are said to 
be antipodal if they are linked by a geodesic, that is if rfoo(C)0 = 2. 
Two sets A^,A+ C dX are antipodal if all pairs (C^O ^ ^ ^+ are 
antipodal. Two measures i^_,z/+ on dX are antipodal when they are 
concentrated on antipodal sets. They are said to be uniformly antipodal 
if their supports are antipodal. 

Note that dX is compact when endowed with the cone topology, so 
that the supports of are compact. If X is a visibility space, then 
antipodality is equivalent to disjointness; in particular two measures are 
then uniformly antipodal if and only if they have disjoint supports. 

Given a complete unit geodesic /i, the asymptotic formula implies that 
Woo(/^-oo, /^+oo) = 2, so that the ends of any complete unit geodesic of 
W^iX) must he antipodal. 

We shall now give a complete answer to Question [1] in the case when X 
is a tree. In particular, we will see that all pairs of uniformly antipodal 
measures, but not all pairs of merely antipodal ones, are the ends of a 
complete geodesic. We shall in fact give two characterizations of ends of 
complete geodesies, but this will need a fair amount of definitions. 

7.2. Flows and antagonism. — From now on, X is assumed to be 
a tree, described as a graph by a couple iy, E) where: V is the set of 
vertices; E is the set of edges, each endowed with one or two endpoints 
in V and a positive length. Since X is assumed to be complete, the edges 
with only one endpoint are exactly those that have infinite length. Since 
X is locally compact, as a graph it is locally finite. It is assumed that 
vertices are incident to 1 or at least 3 edges, so that the combinatorial 
description of X is uniquely determined by its metric structure. We also 
fix a base point xq G X. 
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The following definition is the key to characterize cyclically monotone 
dynamical transport plans. 

Definition 7.2 (antagonism). — We say that two geodesies are an- 
tagonist if there are two distinct points x, y such that one of the geodesies 
goes through x and y in this order, and the other goes through the same 
points in the other order. This amounts to say that there is an edge 
followed by both geodesies, but in opposite directions. 

The point is that for a dynamical transport plan /i, not having pairs 
of antagonist geodesic in its support will prove both a very restrictive 
condition and a characterization of cyclical monotonicity. 

We add to each infinite end a formal endpoint at infinity to unify 
notations. Each edge e has two orientations (xy) and (yx) where x, y 
are its endpoints. The complement in X of the interior of an edge e of 
endpoints x,y has two components Cx{xy) 3 x and Cy{xy) 3 y. An 
orientation, (xy) say, of an edge has a future {xy)+ := Cy{xy) fl dX and 
a past (xy)- := Cx{xy) fl dX. In other words, the future of (xy) is the 
set of all possible endpoints of a complete ray passing through x and y 
in this order. 

Definition 7.3 (Flows). — Assume z/„ and z/+ are antipodal measures 
on dX. Define a signed measure by z/ = z/+ — The flow (defined by 
(z/„,i^+)) through an oriented edge (xy) is 4>{xy) := z/((xy)+). The flow 
gives a natural orientation of edges: an oriented edge is positive if its flow 
is positive, neutral if its flow is zero, and negative otherwise, see figured 
Given a vertex x, let yi, . . . , yk be the neighbors of x such that (xyi) is 
positive, and ^i, . . . , 2;^ be the neighbors of x such that (xzj) is negative. 
Then 4>{xyi) = (t){zjx) is called the flow through x and is denoted 
by 4>{x). If X ^ Xo, then there is a unique edge starting at x along which 
the distance to xq is decreasing. If this edge is a positive one, (xyi^) say, 
then define the specific flow through x as = Xli^io If this 

edge is a negative one, {xzj^^), then let </)°(x) = X^j^jo If ^^i^ edge 

is neutral or if a; = Xq, then let (fp^x) = (pix). 

Note that since is a partition of dX and u has zero 

total mass, we have 

(j)(xy) = -z/((xy)_) = -<p{yx). 

The flow through an edge vanishes if one of the connected component of 
its complement is finite (but this is not a necessary condition). 
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a pendant subtree 



an infinite edge 



Figure 5. The edge (xy) is positive if i^+{{xy)^) > 
since this means that some mass has to move through {xy) 
(in this direction). The edges of the pendant subtree (finite 
complement of a vertex) are ah neutrab 



The next result shows the meaning of flows with regard to our problem. 
Given a dynamical transport plan fi, we denote by fi{xy) the //-measure 
of the set of geodesies that go through an edge (xy) in this orientation, 
by fi{x) the /i-measure of the set of geodesies that pass at x, and by /i°(a;) 
the /i-measure of those that are moreover closest to xq at this time. 



Lemma 7.4- — ///^ is any dynamical transport plan with ends u±, then: 

(1) for all edge {xy) we have [i{xy) ^ max(0(xy), 0), 

(2) for all vertex x we have /i(x) ^ 0(x). 

and each of these inequality is an equality for all {xy), respectively x, if 
and only if /i contains no pair of antagonist geodesies in its support. In 
this case, we moreover have fi^{x) = 0°(x) for all x. 



Note that the set of geodesies that go through an edge {xy) in this 
orientation is open in ^^{X), so that /i contains no pair of antagonist 
geodesies in its support if and only if /i (X> /i-almost no pair of geodesies 
are antagonist. 

Proof. — We only prove the first point, the other ones being similar. 
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Denote by ji{Cy{xy)) the measure of the set of geodesic that he en- 
trirely in Cy{xy). We have 

+ ^i{Cy{xy)) = iy+{{xy)+) = (f){xy) + u_{{xy)+) 

and 

= f^{Cy{xy)) + fi{yx). 
It follows that (j){xy) = fi{xy) — fi{yx) so that fj.{xy) ^ 4>{xy). More- 
over the case of equality fi{xy) = max((/)(x?/), 0) implies that fi{yx) = 
whenever fi{xy) > 0, and we get the conclusion. □ 

From this we deduce in particular a characterization of cyclical mono- 
tonicity. 

Lemma 7.5. — A dynamical transport plan fi is d"^- cyclically monotone 
if and only if ^® ^-almost no pairs of geodesies are antagonist. 

Proof. — Assume that the support of /i contains two antagonist geodesies 
7, (3 and let x, y be points such that = x,'ju = y where u > t and 
= y, = X where w > v. Let r = min(t, v) and s = max(u, w). Then 
(see figure [6]) 

rf(7„ + rf(7„ {3rf < d{^r. 7s)' + d{^, (3sf) 

so that the transport plan {er,es)^fi between /i^ and fig would not be 
cyclically monotone. 




Figure 6. Antagonist pairs of geodesies eontradict the cyclical 
monotonicity: the transport plan corresponding to the solid 
arrow is cheaper than the one corresponding to the dashed 
arrows. 

Assume now that yU ® yU-almost no pairs of geodesies are antagonist. 
Let r : X — 7- M be a function that is continuous, increasing isometric 
on each positive edge and constant on neutral edges. Such a function 
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can be defined locally around any point, and we can design it globally 
since X has no cycle. By Lemma 17.4^ we see that r is isometric when 
restricted to any geodesic in the support of /i (such a geodesic must go 
through positive edges only). Given times r < s, the only | ■ p-cyclically 
monotone transport plan in M from r^fir to r^fig is known to be the 
increasing rearrangement by convexity of the cost. Here r^fis is the r — s 
translate of r#/ir, so that this transport plan has cost (r — s)^. But r is 
1-Lipschitz, so that any transport plan from fir to fig has cost at least 
(r — s)^. This proves the cyclical monotonicity of fi. □ 

Note that here for example, fir and fig need not have finite second 
moment; however fi induces a transport plan with finite cost between 
them, and that peculiarity has therefore no incidence on the proof. 

7.3. Gromov product. — Before we state the main result of this 
section, let us turn to a different point of view that might generalize to 
a wider class of spaces. 

Definition 7.6 (Gromov product). — Given ^_ ^ ^+ in dX, we de- 
note by (^_,^+) C X the locus of a geodesic whose ends are ^_ and 
Then we write 

Do(^,e+) = t^(a:o,(e-,e+)) 
the distance between the base point xq and the geodesic (^_,,^+). Since 
X is a tree, this quantity is equal to what is usually called the Gromov 
product ■ ^+)a;o, see e.g. |BH99] : however the present definition is 
adapted to our needs. Set Dq{^,^) = oo and denote by 7(^_,,^+) the 
parametrized unit complete geodesic whose ends are at ±oo, and 
such that its time realizes Dq: 

c?(xo,7(e-,e+)o) = /^o(^,e+)- 

An important property of the Gromov product is that for any e > 0, 
exp{—eDo) metrizes the cone topology on dX. In particular, -Do('Cn, C) ^ 
oo if and only if — ^ C- Moreover, by compactness of dX, if -Do('Cn.5 Cn) — ^ 
oo then there are increasing indices (rik) such that and Cn^ converge to 
a common point. Another obvious consequence is that Do is continuous, 
and one can even check that it is locally constant outside the diagonal. 
A motivation for using Gromov product here is that it generalizes (with 
some differences though) to a much more general class of spaces, called 
Gromov-hyperbolic, including all CAT(k < 0) Hadamard spaces. One 
could hope to extend part of Theorem I7.1UI below to this class of spaces. 
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Lemma 7.7. — The map 

F -.dX xdX (X) 

is continuous. 

Recall that '^^{X) is endowed with the topology of uniform conver- 
gence on compact subsets. 

Proof. — Since we are dealing with metrizable topologies, we only have 
to check sequential continuity. Assume — )■ ^ and C ^ire two 

converging sequences in dX and set 7 = 7(^, (), 7" = 7(^", C")- 

For any R > 0, for any large enough n, we have that the rays from 
7o to and to ^ coincide on [0, i?]. The same property holds true with 
^'s. The union of (^", 'J-b), {i~b., 1r) and (7^, C") is therefore a geodesic, 
thus must be the locus of 7(^"', C")- Now (^(70, 7q ) is uniformly bounded 
by 2Do{^, C) + 1 for n large enough. It follows that 79 lies on 7, so that 
by definition 70 = 79 and finally 7 and 7" coincide on [—R, R] for n large 
enough. □ 

Since F is a right inverse to (e_oo, e+00), a transport plan 11 G r(z/_, z/+) 
can always be written (e_oo, e+oo)#/^ by taking /i = F^U. We shall 
denote also by Dq the map 7 (-> d{xQ, 7) where 7 is any parametrized or 
unparametrized complete geodesic. 

Lemma 7.8. — A transport plan Uq G r(z/_, z/+) such that 

j -Dl Ho > -00 

is —Dq- cyclically monotone if and only if F^Uq contains no pair of an- 
tagonist geodesies in its support. In this case, Uq is a solution to the 
optimal transport problem 

(6) inf 

U£r{u-,u+) 

Note that the first statement can be recast in the following terms: a 
dynamical transport plan fi with ends i'± and such that / —D^ /i > — 00 
is — DQ-cyclically monotone if and only if it contains no pair of antagonist 
geodesies in its support. 



-Dii^,c)m^,do. 



32 



JEROME BERTRAND & BENOIT KLOECKNER 



Proof. — Consider Ho G r(z/_, and let /i = F^Uq. 

Assume first that there are antagonist geodesies 7, f3 in the support of 
/X. Then permuting 7+00 and /3+oo contradicts the — -DQ-cychcal mono- 
tonicity. 

To prove the other imphcation, assume that supp fi contains no pair 
of antagonist geodesies but that Ho does not achieve the infimum in (|6]). 
This happens notably when Ho is not cyclically monotone. 

Then there is some Hi G r(z/_, z/+) such that 



If F^Ui has couples of antagonist geodesies in its support, then we 
can still improve Hi. Choosing any numbering 62,63,... of the non- 
oriented edges of X, we inductively construct transport plans 112, 1^3, . . . 
in r(i/_,i/_|_) such that -^#11^ has no antagonist geodesies through the 
edges 62,...,6fc in its support, and (— Dq)^!!^ is stochastically domi- 
nated by (— Dq)^!!^-! (for all (CO goii^g through Ck in the negative 
direction, replace ^ by some in the future of 6fc and corresponding to a 
(C',0) ^ suppllfc.i, and replace (C'^O) by (C'lOi there are many choices 
to do but they can be made in a arbitrary manner). Then we get 



where some, or even all of this integrals can be negative infinite. 

We shall use a weak convergence, but — -Dq is not bounded; we therefore 
introduce the functions /r = —min{DQ,T) for all T G N. The transport 
plans Uk also satisfy 



Since dX is compact, so is r(i/_, z/+) and we can extract a subsequence 
of (Ilfc) that weakly converges to some 11, and fi := F^Il has no pair of 
antagonist geodesies in its support. 

The montone convergence theorem implies that for T large enough, 
we have / /rlli < f —Dq Uq, and by weak convergence we get / /rll ^ 
/ /tIIi. Since -D^ ^ fr, we get 






for all T. 
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But then by Lemma [73] we get that = fi^{x) for all x & V, and since 



/ -Dl Ho = Ex- -d^i^o^ it follows that / -D^ Uq and / -D^ ft 



Note that the lemma stays true if —Dq is replaced with any decreasing 
function of Dq, but that we shall need precisely — -Dq loiter. 

Remark 7.9. — In this proof we cannot use Theorem 4.1 of |Vil09] 
since we do not have the suitable lower bounds on the cost. 

7.4. Characterization of ends. — We can now state and prove the 
answer to Question [T] in the case of trees. 

Theorem 7.10. — Assume that X is a tree and let z/_, z/+ be two an- 
tagonist measures on dX . The following are equivalent: 

— there is a complete geodesic in W2{X) with v± as ends; 

— the optimal transport problem (|6]) is finite: 



When these conditions are satisfied, then the optimal transport problem 
([6]) has a minimizer Uq and r^IIo define a geodesic o/^(X) with the 
prescribed ends. 

Moreover, the above conditions are satisfied as soon as u± are uniformly 
antipodal. 

Proof. — First assume that there is a complete geodesic in with 
v± as ends and denote by /i one of its displacement interpolations; by 
Lemma 17.51 the support of /i does not contain any pair of antagonist 
geodesies. From Lemma 17.41 it follows that 



must be equal, a contradiction. 



□ 



ner(i/-,i/+) J 
the specific flow defined by v± satisfies 






since by hypothesis Hq G ^(X). We have also 
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SO that Lemma rTSl implies that Ho := (e_oo, e+oo)#/^ is a solution to 
problem O 

Now consider the case when i/j- are uniformly antipodal. Since the 
supports of v~ and are disjoint, the map -Dq, when restricted to 
supp z/~ X suppz/^, is bounded. Therefore, since it is a continuous map, 
the optimal mass transport problem is well-posed and admits minimisers. 

More generally when the infimum in problem ([6j) is finite, by using 
the regularity of Borel probability measures on dX we can approximate 
v~ and by probability measures whose supports are disjoint sets. 
Then, the previous paragraph gives us a sequence of plans which are 
— Z^Q-cyclically monotone. Since is a continuous map, Prokhorov's 
theorem allows us to extract a converging subsequence whose limit Ho is 
— DQ-cyclically monotone. By the finitness assumption, 

y -Dl{iX) no(rfe,c?C)>-oo 

and Ho is a — Dp-optimal transport plan. 

As soon as a minimizer Ho to ^ exists, by Lemma [7.81 /x := L^IIo is a 
dynamical transport plan that has no antagonist pair of geodesies in its 
support. By Lemma [7.51 /i is cyclically monotone. By its definition 

y c?(x,xo)VoW = J Dl{^X)^oid^.dC) < +00 

so that /i defines a geodesic of W2{X). It has the prescribed ends since 

noer(z/_,z/+). 

We have only left to consider the case when Xlxev (p^{x)d{x, XqY < oo. 
For this, let us construct a suitable complete geodesic by hand. Let r 
be the time function such that t{xq) = 0, as in the proof of Lemma [7.51 
The levels of the time function are finite unions of isolated points and of 
subtrees of X all of whose edges are neutral. Indeed, consider a point a: 
if it lies inside a neutral edge, then all the edge has time r(a). Otherwise, 
let (xy) be the orientation of this edge that is positive: points on [x, a) 
have time lesser than r(a), while points on (a, y] have time greater than 
r(a). If a is a vertex, then similarly one sees that nearby a, only the 
points lying on an incident neutral edge can have time equal to r(a). 
Let f~^{t) be the union of the isolated points of the level T~^{t) and of 
the points that lie (in X) on the boundary of the neutral subtrees of the 
same level. In other words, T~^{t) is the level t of the map induced by r 
on the subforest of X where all neutral edges have been removed. 
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Define now 

where (j){a) = 4>{xy) if a lies inside a positive edge (xy). Note tliat fit is 
a probability measure thanks to the antipodality of z/_ and z/+: without 
it, it would have mass less than 1. It is a good first candidate to be the 
geodesic we are looking for, except the second moment of fit need not 
be finite! To remedy this problem, proceed as follows. First, there is a 
displacement interpolation fi of (fit), which is a probability measure on 
Now, construct a random geodesic 7 as follows : draw 7 with 
law fl, and let 7 be the geodesic that has the same geometric locus and 
the same orientation as 7, and such that 7 is nearest to Xq at time 0. 
The condition ^^.^^ Xq)^ < 00 ensures that the law of 70 has 

finite second moment. □ 

The example shown in Figure [7] shows that antipodality is not sufficient 
for u± to be the ends of a geodesic. 

^ 4 ^ 

^ 2 ^ 

1 



Xo 




d • d • d • d • 

1 11 11 1 1 1 

2 24 48 8 8 8 



Figure 7. The measures (black dots) and (white dots) 
are antipodal. However, the only possible geodesies having 
these measures as ends, depicted by simple arrows, are not in 



PART III 
THE GEODESIC BOUNDARY 



In this part, we adapt to #^(X) the classical construction of the cone 
topology on Hadamard space, see for instance |Bal95j . We introduce this 
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topology on #^(X) = #^(X) U We shall prove in Proposition 

18.11 that the cone topology turns W^iX^ into a first-countable Hausdorff 
space and that the topology induced on coincides with the topol- 

ogy derived from the Wassertein metric. 

In the second part, we shall show in Theorem 19.21 that dW^iX^ is 
homeomorphic to 



endowed with the weak topology. In Corollary 19. 3[ we rewrite the above 
result in terms of Wassertein space over cdX. 

Throughout this part, we will use the following notations. Let (F, d) be 
a geodesic space and y . We set ^y(F) (respectively ^y^i(F)) the set 
of geodesic rays in Y starting at y (respectively the set of unitary geodesic 
rays starting at y\ These sets are closed subsets of ^(F) endowed with 
the topology of uniform convergence on compact subsets. 



The cone topology on W2{X) = >^2(^) UdWi^X) is defined by using 
as a basis the open sets of W2{X) together with 

where x G X is a fixed point, ^ runs over dW2{X), R and e run over 
(0, -|-oo) and fis^^e is the unitary geodesic between 6x and 9 (existence 
and uniqueness follow from Lemma 18. 3p . 

The main properties of the cone topology on y^(X) are gathered to- 
gether in the following proposition. 

Proposition 8.1. — The cone topology on W2{X) is well-defined and 
is independent of the choice of the basepoint 6x- Moreover, endowed with 
this topology, '^2{X) is a first- countable Hausdorff space. By definition, 
the topology induced on W2{X) coincides with the topology derived from 
the Wassertein metric. 

Remark 8.2. — We emphasize that the topology induced on dW2{X) 
by the cone topology coincides with the quotient topology induced by 
the topology of uniform convergence on compact subsets on the set of 
unitary rays in #^(X). 




8. Cone topology on W2{X) 
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Moreover, since (9#2(-^) endowed with the cone topology is first- 
countable, continuity and sequential continuity are equivalent in this 
topological space. 

The scheme of proof is the same as in the nonpositively curved case. 
However, to get the result, we first need to generalize to our setting some 
properties related to nonpositive curvature. 

Lemma 8.3. — Given x E X , the set of unitary rays in starting 
at 5x is in one-to-one correspondence with the set ^i{cdX). Moreover, 
for any ^ G dW2{X), there exists a unique unitary ray starting at 6^ and 
belonging to ^. 

Proof. — Recall that there exists a unique transport plan between a 
Dirac mass and any measure in Since there is a unique geodesic 

between two given points in X, the same property remains true for dy- 
namical transportation plans. Using the previous remarks, we get that 
any fi G ^(^^(X)) such that / s'^{j)^{d'y) < +oo induces a ray start- 
ing at Sx. Moreover, since displacement interpolation always exists (see 
Proposition 14.21) . the set Ms^,i{'p2{X)) is in one-to-one correspondence 
with the unitary dynamical transportation plans starting at 6^, namely 
the measures fi G ^{3?x{X)) such that / s'^{'j)fi{d'y) = 1. Now, since X 
is a Hadamard space, we recall that M{X) is homeomorphic to X x cdX 
(where the left coordinate is the initial location of the ray). Therefore, 
for any x E X, the previous map induces a homeomorphism 

(7) 0. : Mx{X) cdX. 

This gives us a one-to-one correspondence between the set of unitary 
dynamical transportation plans starting at 5x and the set l^i{cdX). 

Given ^ G d'W2{X), consider a unit ray (/if) in 'W2{X) belonging to ^ 
and /ioo e I^i{cdX) its asymptotic measure. We claim that (f>x^#f^oo is 
the unique ray starting at 6^ and belonging to ^. Indeed, (px^^fioo and 
(Ht) have the same asymptotic measure, thus they are asymptotic thanks 
to the asymptotic formula (Theorem 15.21) . The asymptotic formula also 
implies that two asymptotic rays (starting at 6x) have the same asymp- 
totic measure, thus they are equal by the first part of the lemma. □ 

Lemma 8.4- — Let {fit),{(Tt) be two unitary geodesies (possibly rays) 
in W2{X) starting at Sx- Then, for any nonnegative s,t, the comparison 
angle Z^^^^o-t of the triangle A{6x, /is, (Jt) is a nondecreasing function 
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of s and t. Consequently, the map t — y W{lJ't,crt)/t is a nondecreasing 
function as well. 

Proof — We set dm, dg < +00 the length of (/it) and {at) respectively 
and /i, (T the corresponding optimal dynamical plans. Thanks to the 
lifting lemma, we set G r(/i, a) a dynamical plan such that, for given 
s < dm and t < ds, (es,et)#G is an optimal plan. By definition of the 
Wassertein distance, we get, for any s' < s and t' < t, the following 
estimate 

W'(/i.',a,0< f rf^(7(/),7'(O)0(^7,7)- 



Now, the fact that X is nonpositively curved yields 

d\^{s'),y{t')) < f-d\xMs)) + ^^d\x,i{t)) 



s 



t2 



s't' 

- 2—d{x, 7(s))(i(x, 7'(t)) cos Z^(^s)xY{t)- 

where /^.y(s)x'y' (t) is the comparison angle at x (here, we use the fact 
that the initial measure is a Dirac mass). By integrating this inequality 
against 6, we get 



s't' 

+ / c/(x,7(s))d(a;,7'(t)) cos Z^(^)3,y(t) 6(^7, c?7'). 

We conclude by noticing that the inequality above is an equality when 
s = s' and t = t', so we get 

W^(/is',o-t/) < s'^ + - 2s't' cos Zf^^s^^^ 

which is equivalent to the property Z^^,^^^-^, < Z^^s^ac The remaining 
statement follows readily. □ 

Lemma 8.5. — Let (fit) be a unitary geodesic, possibly a ray, starting 
at 6y and 6x 7^ Sy. For any 9 G W2{X) such that 9 7^ 6x, Sy, the compari- 
son angle at 9 satisfies 

cosZs.es, = I d{x,z)d{y,z) cos Z..y9{dz). 

Moreover, given two nonnegative numbers s ^ t, the following inequality 
holds forO<t<T 
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Proof. — For any z E X, the following equality holds 

(f{x, y) = (f{x, z) + (f{y, z) - 2d{x, z)d{y, z) cos ^^zy 

By integrating this inequality against 9, we get the first statement by 
definition of the comparison angle. Let be the unique optimal dynam- 
ical coupling that induces (/if). The first step of the proof is to get the 
equality below: 



cos Z^y^j^^ - 
1 



d{l{t)^l{T))d{'-i{t),x) cos Zx7(t)7(T) Kdl) 



For 7 G supp /i, the following equality holds 

d\^{T),x) = 

d\^{t),x) + d\^{t),^{T)) - 2rf(7(t), x)d(7(t), 7(T)) cos 2,^(i)^(T). 

By integrating this equality against yU, we get ([8]). Now, using that X is 
nonpositively curved, we have ^x-y(t)'y(T) + ^x-t(t)y > ^r. This gives 

cos Z^y^j5^ 

< 



< 



w(MT,MOw(Mt,^.) / d{-i{t),-i{T))d{-i{t),x) cos Z,^(t)j, /i(c?7) 
(T-t)wU,3.) / ^d{-f{t), y)d{-fit), x) cos Z^^(t)j, fx{dj) 
^ wM')w{f^tA) yMlit),x) cos Z,,^(t)j, /i(d7) 

< - COS Z^oM,5^ 

where the last inequality follows from the first statement and the result 
is proved. □ 

As a consequence, we get the following result. 

Proposition 8.6. — Given e > 0, a > 0, and R > 0, there exists 
a constant T = T{e, a,R) > such that the followings holds: for any 
x,y E X such that d{x,y) = a and a unitary geodesic (possibly a ray) 
(fit) of length greater than T and starting at 6y, if (af) is the unitary 
geodesic from 6x to fig then 

WK,a^') <e 
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for any s' > s > T . 

In particular, if [fit) is a ray and s goes to infinity, {cr^)s>o converges 
uniformly on compact subsets to the unitary ray fis^^^ where C, is the 
asymptote class of (fit)- 




Figure 8. Uniform convergence of on compact subsets. 

Proof — Under these assumptions, the comparison angle Z^^^^^^ is ar- 
bitrary small provided s is sufficiently large. Consequently, thanks to 
Lemma Z^^,^^^^ is close to vr; therefore the comparison angle Z^^^^,^^, 
is small. This gives the ffist part of the result since 

thanks to Lemma 18. 4[ Using Lemma 18.41 again, it only remains to prove 
the pointwise convergence of {a'^)s>o to fis^,^- Thanks to the asymptotic 
formula, there exists C > such that 

(9) W((/i5.,^)t,/it) <C 

for any nonnegative number t. Finally, we conclude by using s' sufficiently 
large and the bound 

where the same reasoning as above and ^ show that the first term on 
the right-hand side is small provided s' is large. □ 

Now, we can prove that the topology above is well-defined and does not 
depend on the choice of the base point 5^. This the content of the lemma 
below. 



GEOMETRIC WASSERSTEIN SPACES 



41 



Lemma 8.7. — Given two positive numbers R,e and y G X , ^ E 
U{x,r],R,e) ndW2{X), there exists S,e' >0 such that 

U{y,^,S,e')cU{x,r],R,e). 

Proof. — We set a = e - W{{fi5^,n)R, > 0. Let 9 e U {y, ^, S, e') 

and 9 (respectively S) be the unitary geodesic fiSy,d (respectively the 
unitary ray iJ^Sy,^). We have 

W((/i5„e)ii, (/i5,,,;)ii) < 

W((/i5„0)ij, (/^5,,es)i?) + W((/i5„es)/?, {.^^8...^s)R) 

+ W((/i5,,Hs)iJ, {^^S.,^R) + W((/i5,,€)i?, (/i5,,,;)i?) 

The first and the third term on the right-hand side are smaller than a/S 
for large S thanks to Proposition 18.61 while the second term is smaller 
than a/S for large S and small e' thanks to lemma 18^ □ 

9. The boundary of W-2{X) viewed as a set of measures 

To state the main result of this part, we first need to introduce a 
definition. 

Definition 9.1. — We set 

Am: ^5^,i(^2(X)) ^i(c9X) 

(/it) I > Atoo 

the map that send a unitary ray starting at 5x to its asymptotic measure. 

The main result of this part is the following theorem. 

Theorem 9.2. — The map Am : ^i^^i{W2{X)) — > ^i{cdX) induces 
a homeomorphism from dW2{X) onto J^i{cdX). 

Note that a straightforward consequence of the result above is 

Corollary 9.3. — Let d he a metric on dX that induces the cone topol- 
ogy on dX and dc the cone metric induced by d on cdX (see ^ for a 
definition). Let us denote by W2{cdX) the quadratic Wassertein space 
over the Polish space {cdX,dc). Then, dl^2{X) is homeomorphic to the 
subset of probability measures with unitary speed in W2{cdX). 
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Remark 9.4- — In particular we get the more symmetric result that 
cd'W2{X) is homeomorphic to W2{cdX). 

The rest of this part is devoted to the proof of the theorem above. 
Recall that we have proved in Lemma 18.31 that both Am and the map 
Am : dW^iX) — >■ !^x{cdX) it induces are bijective. 

The proof of Theorem l9.2l is in two steps. First, we prove that the map 
Am is a homeomorphism. Then, we use this fact to prove that Am is a 
homeomorphism as well. 

We start the proof with a definition. 

Definition 9.5. — Let x E X. We denote by 

ODT, = {fie^{M,{X))-js\^)fi{dj) = l} 

the set of unitary dynamical transport plans endowed with the weak 
topology. We also set 

(ei)>o#: ODT, — > ^5.,i(#^2W) 
I — > (fit) 

and 

eoo#: ODT, ^,{cdX) 
/i I — y 

where (p, is defined in ([7]). 

Thanks to Lemma 18.31 we have the following commutative diagram 
where all the maps are one-to-one. 

ODT, "^-^ » ^i{cdX) 

Am ^-"^ 

(et)t>o# I 

M5^,,{W2{X)) 

We first prove that 

Lemma 9.6. — The map eoo# is a homeomorphism onto ^i{cdX). 

Proof. — The map (j), : J^x{X) — y cdX is a homeomorphism. There- 
fore it induces a homeomorphism between ODT, and ^i{cdX) when 
endowed with the weak topology. □ 

Lemma 9.7. — The map {et)>o^ is a continuous map. 
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Proof. — Since the spaces we consider are metrizable, we just have to 
prove the sequential continuity. Consequently, we are given a sequence 
(A*n)neN such that /i„ ^ /i in ODT^. Now, since Ct : ^xi.^) — >■ is a 
continuous map, we get that et^fin ^t#f^ in ^{X). By definition of 
ODTx, we have 



Since J 3^(7) yU„((i7) = J 7(1)) /i„((i7), the equality above implies 
the convergence of the second moment. Namely, we have 

(f{x,-f{t))nn{d-f)=t'^j (f{x,'j{l))nn{d^) = J rf^(x,7(t)) /i((i7). 

This implies the convergence of et^fin to Cj^/i with respect to the Was- 
sertein distance (see for instance |Vil09] ). Thus, we have proved the 
pointwise convergence of rays. Now, since t — > W{et:ff,fin, is non- 

decreasing as proved in Lemma [8.41 we get the result. □ 

We end the first part of the proof with the following lemma. 

Lemma 9.8. — The map {et)>o^ is a homeomorphism. 

Proof. — Since the topology of both and ODT^ is induced 

by a metric, it is sufficient to prove that (et)>o^ is a proper map. More- 
over, we just have to prove sequential compactness. We set K a compact 
subset of Ms^,i{W2{X)). Let (/U„)„gN ^ {et)>Q^{K). We first notice that 
{lJLn)nm is tight. Indeed, by assumption on K, the sequence {ei^iin)n& 
is tight in !^{X). Therefore, by arguing as in the end of the proof of 
Proposition 14. 2^ we obtain the claim. Consequently, since Mx{.X) is a 
Polish space, we can apply Prokhorov's theorem to get a converging sub- 
sequence {fink)km to ]i. It remains to prove that Jl G ODT^, namely 
that J s"^ {'-f) Jl{d'y) = 1. Since K is compact, we can also assume with- 
out loss of generality that {et^iin^) — > {fit) in 3?Sj:,i{'^2{X)). Moreover, 
since et^fin^ ^ for any t, we get (et)>o^(/x) = {fit)- Therefore, 

ei^fin^ — >■ ei^Ji in #^(X). This implies the convergence of the second 
moment / d'^{x,'j{l)) finkid'j) = J s'^^'j) fin^^d'j) = 1 (see for instance 
|Vil09] ) and the result is proved. □ 
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We are now in position to prove Theorem 19.21 We set pa^^j ^^e canoni- 
cal projection on 9#2(^)- We have the following commutative diagram. 

We have seen at the beginning of the proof that all the maps above 
are one-to-one. To conclude, it remains to prove that is a con- 

tinuous map. Since dW2{X) is first-countable (see Remark 18. 2p . it is 
sufficient to prove sequential continuity. To this aim, let — > C, in 
dW2{X) and {et#fin), (/it) e e^5^,i(#^(X)) such that Pa^2((et#/in)) = 
and PdW2{{f^t)) = ^- Recall that under these assumptions, the map 
t — >■ W{et#fin, fJ't) is nondecreasing (see Lemma . thus we just have 
to show the pointwise convergence of (et^fin)- This pointwise conver- 
gence follows readily from the definition of the cone topology on W2{X). 




PART IV 
ISOMETRIC RIGIDITY 



10. Statement of the result and reduction to a Radon 

transform 

In this part we study the isometry group of W2{X). We have a natural 
morphism 

#:IsomX -> lsomW2{X) 
^ <^#- 

An isometry $ of #2 (A) is said to be trivial if it is in the image of 
to preserve shape if for all /i G W2{X) there is an isometry ip (depending 
on /i) of A such that $(/i) = f#lJ', and to be exotic otherwise. When all 
isometries of W2{X) are trivial, that is when 7^ is an isomorphism, we 
say that 'W2{X) is isometrically rigid. 

The main result of this part is the following. 
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Theorem 10.1 (isometric rigidity). — Let X he a geodesically com- 
plete, negatively curved Hadamard space. If 

(1) X is a tree, or 

(2) X is a Riemannian manifold, or more generally 

(3) there is a dense subset of X all of whose points admit a neighbor- 
hood isometric to an open subset of a Riemannian manifold of dimension 
at least 2, 

then #2(^) is isometrically rigid. 

By "Negatively curved" we mean that X is not a line and the CAT(O) 
inequality is strict except for triangles all of whose vertices are aligned. 
Note that when X is Euclidean, is not isometrically rigid as proved 

in |Klo08] . The Ip^q building, as well as many polyhedral complexes 
whose faces are endowed with hyperbolic metrics, provide examples of 
isometric rigidity covered by the third case of the Theorem. 

The strategy of proof is as follows: we show first that an isometry must 
map Dirac masses to Dirac masses. The method is similar to that used in 
the Euclidean case, and is valid for geodesically complete locally compact 
spaces without curvature assumption. It already yields Corollary 110. 5[ 
answering positively in a wide setting the following question: if #2(^) 
and W2{Z) are isometric, can one conclude that Y and Z are? 

The second step is to prove that an isometry that fixes all Dirac masses 
must map a measure supported on a geodesic to a measure supported 
on the same geodesic. Here we use that X has negative curvature: this 
property is known to be false in the Euclidean case. 

Then as in |Klo08] we are reduced to prove the injectivity of a spe- 
cific Radon transform. The point is that this injectivity is known for 
MH" and is easy to prove for trees, but seems not to be known for other 
spaces. We give a partial but sufficient injectivity argument for mani- 
folds. This argument can be generalized to spaces with a dense set of 
higher- dimensional manifold-type points. 

10.1. Characterization of Dirac masses. — The characterization 
of Dirac masses follows from two lemmas, and can be carried out without 
any assumption on the curvature. Note that when we ask a space to be 
geodesically complete, we of course imply that it is geodesic. 
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Lemma 10.2. — Let Y he a locally compact, geodesically complete Pol- 
ish space. Any geodesic segment in #2(^) issued from a Dirac mass can 
be extended to a complete geodesic ray. 

Proof. — Let us recall the measurable selection theorem (see for example 
|Del75] . Corollary of Theorem 17): any surjective measurable map be- 
tween Polish spaces admits a measurable right inverse provided its fibers 
are compact. Consider the restriction map 

: M{Y) ^ ^O'^(F) ; 

the fiber of a geodesic segment 7 is closed in the set e^^Q,s(F) of geodesic 
rays of speed s = 3(7) starting at 70. This set is compact by Arzela- 
Ascoli theorem (equicontinuity follows from the speed being fixed, while 
the pointwise relative compactness is a consequence of Y being locally 
compact and geodesic, hence proper). 

Moreover the geodesic completeness implies the surjectivity of p^. 
There is therefore a measurable right inverse q of p'^: it maps a geodesic 
segment 7 to a complete ray whose restriction to [0,T] is 7. 

Let {fit)t£io,T] be a geodesic segment of #2(^), of speed s, with fiQ = Sx 
and let = q^{^T)'- it is a probability measure on ^{Y) whose restriction 
p'^ifi) is the displacement interpolation of (fit). For all t > we therefore 
denote by fit the measure {et)#fi on X. 

It is easy to see that it defines a geodesic ray: since there is only one 
transport plan from a Dirac mass to any fixed measure, W^(A*o!/^t) = 
/ s{'y)H'^ ^{d'j) = sH"^. Moreover the transport plan from fit to fit' de- 
duced from ft gives W(/it,/it') ^ s\t — t'\. But the triangular inequality 
applied to fiQ, fit and fit' implies W(yUt, yUf) ^ s\t—t'\ and we are done. □ 

Lemma 10.3. — If Y is a geodesic Polish space, given fio € #2(^) 'not 
a Dirac mass and y G supp fiQ, the geodesic segment from fiQ to fii := 5y 
cannot he geodesically extended for any time t> 1. 

Note that at least in the branching case the assumption y G supp /iq is 
needed. 

Proof. — Assume that there is a geodesic segment {fit)t&[Q,i+e\ and let 
fi be an optimal dynamical transport plan from /iq to fii+£. Since y G 
supp /io, supp fi contains a geodesic segment of X that is at y at times 
and 1, and must therefore be constant. 

Since fiQ is not a Dirac mass, there is a y' 7^ y in supp fiQ. Let 7 be a 
geodesic segment in supp fi such that 70 = y' . Then 71 = y lies between 
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y' and y" := 71+^ on 7, and the transport plan 11 = ^^'^^^ contains in its 
support the couples {y",y') and {y,y). 

But then cyclical monotonicity shows that 11 is not optimal: it costs 
less to move y" to y and y to y' by convexity of the cost (see figurelH]). The 
dynamical transport fi cannot be optimal either, a contradiction. □ 




Figure 9 . The transport shown with continuous arrows is less 
effective than the transport given by the dashed arrows. 

We can now easily draw the consequences of these lemmas. 

Proposition 10.4- — Let Y, Z be geodesic Polish spaces and assume 
that Y is geodesically complete and locally compact. Any isometry from 
W2(Y) to W2{Z) must map all Dirac masses to Dirac masses. 

Except in the next corollary, we shall use this result with Y = Z = X . 
Note that we will need to ask that X be geodesically complete in addition 
to the Hadamard hypothesis. 

Proof. — Denote by ip an isometry — )■ WiiZ") and consider any 
X ^ Y . If i^i^x) were not a Dirac mass, there would exist a geodesic 
segment (/if )fe[o,i] from ^{bx) to a Dirac mass (at a point y G supp V5(5x)) 
that cannot be extended for times t > 1. 

But ip~^{fj,t) gives a geodesic segment issued from 6x, that can therefore 
be extended. This is a contradiction since (p is an isometry. □ 

Corollary 10.5. — Let Y, Z be geodesic Polish spaces and assume that 
Y is geodesically complete and locally compact. Then #2(^) is isometric 
to W2{Z) if and only ifY is isometric to Z. 

Proof. — Let ip be an isometry W2{Y) — #^(Z). Then ip maps Dirac 
masses to Dirac masses, and since the set of Dirac masses of a space is 
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canonically isometric to the space, 93 induces an isometry Y ^ Z. the 



Note that we do not know whether this result holds for general metric 
spaces. Also, we do not know if there is a space Y and an isometry of 
that maps some Dirac mass to a measure that is not a Dirac mass. 

10.2. Measures supported on a geodesic. — The characterization 
of measures supported on a geodesic relies on the following argument: 
when dilated from a point of the geodesic, such a measure has Euclidean 
expansion. 

Lemma 10.6. — Assume that X is negatively curved, and let 'j be a 
maximal geodesic of X , fi be in #^(X). Given a point x G 7, denote by 
(x* ■ /i)tg[o,i] the geodesic segment from 6^ to fi. If y is any point of X, 
{x + y)/2 denotes the midpoint of x and y. 

The measure fi is supported on 7 if and only if for all x,g & 'j 



Proof. — The "only if" part is obvious since the transport problem on 
a convex subset of X only involves the induced metric on this subset, 
which here is isometric to an interval. 

To prove the "if" part, first notice that the CAT(O) inequality in a 
triangle {x,y,g) yields the Thales inequality d{{x + y)/2, (x + g)/2) ^ 
^d{y,g), so that by direct integration 



Assume now that /i is not supported on 7 and let y G supp /i \ 7, and 
x,g & 'J such that x, y, g are not aligned. Since 7 is maximal, this is 
possible, for example by taking d{x, g) > d{x, y) in the branching case 
(any choice oi x ^ g would do in the non-branching case). Since X is 
negatively curved, we get that d{{x,y)/2, {x + g)/2) < \d{y,g) so that 
ffTOj) is strict (see Figure [JO]). □ 

Corollary 10.7. — If X is a negatively curved Hadamard space and 
is a maximal geodesic of X , any isometry of W2{X) that fixes Dirac 
masses must preserve the subset #2(7) of measures supported on 7, and 
therefore induces an isometry on this set. 



converse implication is obvious. 



□ 





W(X2 ■ /i, (5(^.+g)/2) ^ -W(/i, 5g). 
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Figure 10. In negative curvature, midpoinds are closer than 
they would be in Euclidean space. 



Proof. — Let (p be an isometry of W2{X) that fixes Dirac masses. Let 
7 be a maximal geodesic of X and /i G W2{X) be supported on 7. 
If (p{fi) where not supported in 7, there would exist x, (7 G 7 such 
that W(v5(/x)i/2, 5(x+g)/2) < I W(v5(/i), But ip is an isometry and 

P'-^iSx) = Sx, P'-^iSg) = 6g, (p~^{6(x+g)/2) = S(x+g)/2 SO that alsO V5(/i)l/2 = 

As a consequence, ffTOl) would be a strict inequality too, a con- 
tradiction. □ 

10.3. Isometry induced on a geodesic. — We want to deduce from 
the previous section that an isometry that fixes all Dirac masses must also 
fix every geodesically-supported measure. To this end, we have to show 
that the isometry induced on the measures supported in a given geodesic 
is trivial, in other terms to rule out the other possibilities exhibited in 
|Klo08] at which it is recommanded to take a look before reading the 
proof below. 

Proposition 10.8. — Assume that X is a geodesically complete, neg- 
atively curved Hadamard space and let ip he an isometry o/^(X) that 
fixes all Dirac masses. For all complete geodesies 7 of X , the isometry 
induced by p on #2(7) is the identity. 

We only adress the geodesically complete case for simplicity, but the 
same probably holds in more generality. 

Proof. — Let x be a point not lying on 7. Such a points exists since X 
is negatively curved, hence not a line. 

First assume that p induces on #2(7) an exotic isometry. Let y, z be 
two points of 7 and define /xq = + \^z- Then := p"'{fio) has the 
form rrinSy^ + (1 - mn)S^^ where m„ -» 0, y„ 00, z„ -> + z)/2 =: g. 

Let now 7' be a complete geodesic that contains y' = {x + y)/2 and 
z' = {x + z)/2 (see Figure [TT]). Since the midpoint of /io and 6^ is 
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supported on 7', so is the midpoint of /i„ and 5^. This means that 
{x + Un)/'^ and (x + Zn)/^ he on 7', and the former goes to infinity. 
This shows that also induces an exotic isometry on #2(7')' thus that 
{x + Zn)/2 goes to the midpoint g' of y' and z' . But we already know 
that (x + 2n)/2 tends to the midpoint of x and which must therefore 
be g'. 



X 




Figure 11. That an isometry acts exotically on a geodesic 7 
would imply that there is another geodesic 7' such that taking 
midpoints with x is an affine map. Then 7 and 7' would be 
parallel, therefore ound a flat strip. 



Since this holds for all choices of y and 2, we see that the map y — )■ 
{x + y)/2 maps affinely 7 to 7'. But the geodesic segment [x7t] converges 
when t — 7- ±00 to geodesic segments asymptotic to 7 and —7 respectively. 
It follows that 7' is parallel to 7, so that they must bound a flat strip (see 
jBal95] ). But this is forbidden by the negative curvature assumption. 

A similar argument can be worked out in the case when induces an 
involution: given y, z ^ and their midpoint g, one can find measures 
fin supported on y and ?/„, where yn ^ 9 and yU„ has more mass on yn 
than on y, such that </?(/in) is supported on z and a point Zn of 7, with 
more mass on Zn- It follows that the midpoints y\y'ni z' with x are on a 
line, and we get the same contradiction as before. 

The classification of isometries of ^(M) shows that if (/? fixes Dirac 
masses and is neither an exotic isometry nor an involution, then it is the 
identity. □ 

Now we are able to link the isometric rigidity of W^i/Oj to the injectiv- 
ity of a Radon transform. The following definition relies on the following 
observation: since a geodesic is convex and X is Hadamard, given a point 
y and a geodesic 7 there is a unique point v-^iy) G 7 closest to called 
the projection of ?/ to 7. 
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Definition 10.9. — When X is geodesically complete, we define the 
perpendicular Radon transform i^/i of a measure /i G W2{X) as the 
following map defined over complete geodesies 7 of X: 

In other words, this Radon transform recalls all the projections of a 
measure on geodesies. 

The following is now a direct consequence of Proposition 110. 81 

Corollary 10.10. — Assume that X is geodesically complete and neg- 
atively curved. If there is a dense subset A C W2{X) such that for all 
fi G W2{X) and all u & A, we have 

then #2(X) is isometrically rigid. 

Proof. — Let (f be an isometry of W2{X). Up to composing with an 
element of the image of ^, we can assume that (p acts trivially on Dirac 
masses. Then it acts trivially on geodesically supported measures. For 
all u E A, since {p-y)#i^ is the measure supported on 7 closest to u, one 
has ^(p^i^) = ^i', hence ip^u) = p. We just proved that (/? acts trivially 
on a dense set, so that it must be the identity. □ 

11. Injectivity of the Radon transform 

The proof of Theorem 110.11 shall be complete as soon as we get the 
injectivity required in Corollary 110.101 Note that in the case of the real 
hyperbolic space MH", we could use the usual Radon transform on the 
set of compactly supported measures with smooth density to get it (see 
|Hel99j ). 

11.1. The case of manifolds and their siblings. — Let us first give 
an argument that does the job for all manifolds. We shall give a more 
general, but somewhat more involved argument afterward, and the next 
subsection shall deal with the case of trees. 

Proposition 11.1. — Assume that X is a Hadamard manifold. Let A 
be the set of finitely supported measures. For all fi G W2{X) and all 
V e A, if ^% n = then /i = z/. 
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Note that A is dense in #^(X) so that this proposition ends the proof 
of Theorem 110.1^ case [21 

Proof. — Write u = Yl'^^i^xi where Y^frii = 1. Note that since X is a 
negatively curved manifold, it has dimension at least 2. 

First, we prove under the assumption 1% [i = that /i must be sup- 
ported on the Xj. Let x be any other point, and consider a geodesic 7 
such that 7(0) = X and 7(0) is not orthogonal to any of the geodesies 
[xxi). Then for all 'p-^ixi) 7^ x and there is an e > such that the 
neighborhood of size e around x on 7 does not contain any of these pro- 
jections. It follows that ^z/(7) is supported outside this neighborhood, 
and so does ^/i(7). But the projection on 7 is 1-Lipschitz, so that \i 
must be supported outside the e neighborhood of x in X. In particular, 
X ^ supp /i. 

Now, if 7 is a geodesic containing Xj, then ^1^(7) is finitely supported 
with a mass at least at Xj. For a generic 7, the mass at Xj is exactly 
mj. It follows immediately, since \i is supported on the Xj, that its mass 
at Xi is mi. □ 

The above proof mainly uses the fact that given a point x and a finite 
number of other points Xj, there is a geodesic 7 9 x such that P'y(xj) 7^ x 
for all i. It follows that the proof can be adjusted to get the case [3] of 
Theorem llO.il To this end, let us define the regular set of X as the set of 
points having a neighborhood isometric to an open set in a Riemannian 
manifold of dimension at least 2. 

Proposition 11.2. — Assume that X is a Hadamard, geodesically com- 
plete space. Let A be the set of measures supported on finite set of points 
all located in the regular set of X. For all fi G W2{X) and all v G A, if 
^/i = ^z/ then /i = z/. 

Proof. — Write u = Yl^i^Xi where ^rrij = 1. Consider one of the x,, 
and let := fi{{xi}). We get m'i < rrii from ^ = v. 

Since by hypothesis, Xj has a neighborhood Ui isometric to an open 
set in a Riemannian manifold of dimension at least 2, as before there is 
a geodesic 7 9 Xi such that p-y^Xj) 7^ Xj for all j 7^ i. It follows that 
for some e > 0, the measure M fJ'il) = ^ 1^(7) is concentrated outside 
B{xi,e) \ Xi. By taking e small enough so that B{xi,e) is contained in 
Ui, it follows also that /i is concentrated outside B{xi,e) \ Xj. 

Let 7-*- denote the set of points x such that the geodesic segment [xxj] 
is orthogonal to 7 at Xj. Then 7-*- fl B{xi,e) is the exponential image of 
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the orthogonal of the tangent hne of 7 at x^. By definition of the Radon 
transform, /i(7"'" \ Xj) = — and for all x G supp fi \ 7-*", we have 
p^{x) ^ B{xi,e). 

Choose a second geodesic 72 3 Xt, close enough to 7 to ensure that 
for all X G supp/i \ 7"*", we still have P72(^) ^ B{xi,e) (up to schrinking 
e a bit if necessary). We can moreover assume that 72 enjoys the same 
properties we asked to 7, so that 

Vx e supp /i \ 7"*" n 7^, p^2{x) ^ B{xi,e). 

Let n be the dimension of f/j. We can construct inductively a fam- 
ily 71, .. . ,7„, of geodesies chosen as above, and such that their velocity 
vectors at Xi span its tangent space. For all point x in 

supp/i\ na7a> 

we get p.y„{x) ^ B{xi,e). But Ciala = {^i}^ and considering ^/i(7„) = 
^ J^iln) we get m'i = rrii. Since /i is a probability measure and '^nii = 1, 
we deduce = v. □ 

Now if X fits in case [3] of Theorem 110.11 then the regular set of X is 
dense, thus A is dense in #^(X). We can therefore conclude by Corollary 

mm 

11.2. The case of trees. — Let X be a locally finite tree that is not 
a line. Recall the notations already used in Section [61 We describe X by 
a couple (V, E) where V is the set of vertices; E is the set of edges, each 
endowed with one or two endpoints in V and a length. 

For all X & V , let k{x) be the valency of x, that is the number of edges 
incident to x. We assume that no vertex has valency 2, since otherwise 
we could describe the same metric space by a simpler graph. Note also 
that edges with only one endpoint have infinite length since X is assumed 
to be complete. 

In this setting, X is geodesically complete if and only if it has no leaf 
(vertex of valency 1). Assume this, and let 7 be a complete geodesic. If 
X is a vertex lying on 7 and Ci, . . . , Ck{x) are the connected components 
of X \ {x}, let -L^ (7) be the union of x and of the Ci not meeting 7 (see 
Figure [T2|) . It inherits a tree structure from X. In fact, -L^(7) depends 
only upon the two edges e, / of 7 that are incident to x. We therefore 
let ±^(e/) = ±^(7). 

The levels of are called the perpendiculars of 7, they also are the 
bissectors of its points. They are exactly: 
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■ 7 



Figure 12. A perpendicular to a geodesic. 



— the sets -L^(7) where x is a vertex of 7. 

— the sets {x} where a: is a point interior to an edge of 7. 

A measure fi G W2{X) can be decomposed into a part supported out- 
side vertices, which is obviously determined by the projections of fi on 
the various geodesies of X, and an atomic part supported on vertices. 
Therefore, we are reduced to study the perpendicular Radon transform 
reformulated as follows for functions defined on V instead of measures 
on X. 

Definition 11.3 (combinatorial Radon transform) 

A flag of X is defined as a triple (x, e/) where x is a vertex, e ^ f are 
edges incident to x and ef denotes an unordered pair. Let us denote the 
set of flags by F; we write a: G e/ to say that (x, ef) G F. 

Given a summable function h defined on the vertices of X, we define 
its combinatorial perpendicular Radon transform as the map 

^h: F ^ R 

{x,ef) ^ J2 ^(y) 

?/6±-(e/) 

where the sum is on vertices of _L^(e/). 

It seems that this Radon transform has not been studied before, con- 
trary to the transforms defined using geodesies |BCTCP91j . horocycles 
|Car73|. IBFP89j and circles |CTGP03j . 

Theorem 11. 4 (Inversion formula). — Two maps h,l : V ^ R such 
that Yl^ ~ ^'^^ M h = are equal. More precisely, we can recover 
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h from ^ h by the following inversion formula: 

Hx) = E ^ h{x, ef) - My) 

where the first sum is over the set of pairs of edges incident to x. 
Proof. — The formula relies on a simple double counting argument: 
Y,^%h{x,ef) = J2 

vev 

where nx{y) is the number of flags (x, ef) such that y G l-^{ef). \iy ^ x, 
let Cy be the edge incident to x starting the geodesic segment from x to 
y. Then y G ±^(e/) if and only if e, / 7^ Cy. Therefore, nx{y) = (^''^2~^)' 
But n^{x) = (^f^), so that 

E ^ = ^ ^^^^ + ^^^""^ " 

□ 
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